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Abstract 


The  infrared- laser  excitation  of  the  internal  vibrational  mode  of  a 
diatomic  molecule  adsorbed  on  a  metal  surface  is  analyzed  theoretically. 

This  vibrational  energy  is  damped  into  the  metal  by  electron-hole 
excitations.  Simple  expressions  for  the  populations  of  the  vibrational 
levels,  the  mean  number  of  vibrational  quanta  and  the  rate  of  energy 
transfer  between  the  infrared  laser  and  the  metal  surface  at  the  steady 
state  are  derived.  An  equation  of  evolution  can  readily  be  solved 
numerically  to  determine  the  time  necessary  to  reach  this  steady  state.  The 
criteria  of  applicability  of  the  Markov  approximation  (which  leads  to  the 
golden  rule)  is  clearly  established,  where  it  is  seen  that  this 
approximation  may  not  be  used  to  compute  the  evolution  of  the  populations  of 
the  vibrational  levels.  The  random-phase  approximation  is  shown  to  give  the 
correct  kinetic  equation  for  the  populations  of  the  vibrational  levels.  The 
excitation  of  carbon  monoxide  adsorbed  on  a  copper  surface  is  analyzed 
quantitatively. 
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I.  Introduction 


The  excitation  of  a  molecular  vibration  can  significantly  contribute  to 
overcoming  the  activation  barrier  of  a  chemical  reaction.  This  fact  has  led 
to  many  attempts  to  modify  the  course  of  surface  processes  by  resonant 
excitation  of  vibrational  modes  by  means  of  infrared  laser  radiation  [1]. 

Due  to  their  importance  in  catalysis,  laser -stimulated  reactions  at  metal 
surfaces  have  been  investigated  [2-6].  For  example,  it  has  been  shown  that 
focused  infrared  laser  radiation  can  induce  resonant  ionization  of  adspecies 
[2]  or  gaseous  species  above  a  surface  [3],  where  the  fragments  from  the 
latter  process  can  then  undergo  surface  reactions  which  are  not  possible 
with  the  original  molecules.  Even  when  the  laser  is  not  focused,  the 
resonant  excitation  of  a  gaseous  species  can  decrease  the  likelihood  of 
physisorption  [A]  or  increase  the  tendency  for  dissociative  chemisorption 
[5].  When  laser  radiation  impinges  on  a  gas-solid  interface,  it  can  excite 
the  substrate,  adsorbate  and/or  species  in  the  gas  phase.  Chuang  [7]  has 
demonstrated  that  resonant  excitation  of  the  adsorbate  can  stimulate  etching 
of  a  silicon  surface  by  SF^.  While  this  mechanism  of  adsorbate  resonant 
excitation  has  been  proposed  for  experiments  involving  metal  surfaces  [6], 
it  is  perhaps  not  as  obvious  here  since  the  vibrational  energy  of  a  metal 
adsorbate  quickly  damps  into  the  metal,  as  observed  by  infrared  reflection 
absorption  spectroscopy  [8,9].  This  leads  to  two  questions:  (1)  Can  an 
infrared  laser  beam  maintain  a  high  surface  concentration  of  vibrationally- 
excited  molecules,  and  (2)  what  is  the  rate  of  energy  transfer  from  the 
laser  to  the  metal  through  the  adsorbed  molecules?  Such  energy  transfer 
leads  to  resonant  heating  of  the  substrate,  which  can  then  "thermally" 
enhance  a  surface  reaction  or  the  desorption  of  the  adsorbate  before 
reaction  occurs.  These  questions  have  motivated  us  to  carry  out  a 


theoretical  study  of  the  simple  case  of  the  resonant  interaction  of  infrared 

laser  radiation  with  a  diatomic  molecule  adsorbed  on  a  metal  surface. 

The  spectral  infrared  absorption  width  of  a  diatomic  molecule  in  the 
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gas  phase  is  of  the  order  of  10  cm  .  When  the  molecule  is  adsorbed  on  a 

metal,  this  width  may  reach  many  cm  *  [8).  This  broadening  was  successfully 

explained  by  Persson  in  terms  of  electron  transfer  between  the  adsorbate  and 

the  substrate  [10],  and  this  mechanism  has  been  confirmed  by  means  of 

cluster  calculations  [11].  The  lifetime  of  an  excited  vibrational  state  was 

-12 

estimated  to  be  of  the  order  of  10  s.  This  is  much  shorter  than  the 

lifetime  due  to  other  mechanisms  such  as  phonon-phonon  energy  transfer 

_Q  .](] 

(-10  s)  [12-17]  or  dipole-hole-electron  interactions  (-10  s)  [18].  We 
shall  thus  assume  electron  transfer  to  be  the  only  mechanism  responsible  for 
energy  tranfer  between  the  substrate  and  the  adsorbate. 

Theoretical  studies  already  presented  in  the  literature  related  to  the 
interaction  of  an  infrared  laser  beam  with  adsorbed  molecules  are  focused  on 
the  phonon-phonon  energy  transfer  between  the  substrate  and  the  adsorbate. 
The  computational  techniques  used  are  the  "golden  rule"  [13,14]  and  the 
transformation  to  dressed  states  [15-17].  The  Zwanzig  projector  technique 
[19]  within  the  Born,  the  random-phase  and  the  Markov  approximations  has 
been  used  to  compute  molecular  desorption  rates  [12].  We  shall  use  the 
Zwanzig  projector  technique  within  the  Born  approximation  only,  and  shall 
obtain  simple  expressions  describing  the  kinetic  and  the  steady  state  of  our 
system.  Furthermore,  our  results  will  provide  criteria  for  the 
applicability  of  all  the  approximations  listed  above. 

In  Section  II  we  present  our  model  and  write  down  the  Hamiltonian.  In 
Section  III  we  derive  the  kinetic  equation  for  the  density  operator 
characterizing  the  vibrational  state  of  the  adsorbed  molecule.  In  Section 


IV  we  analyze  the  memory  kernels,  where  the  validity  of  the  Markov 
approximation  is  questioned.  In  Section  V  we  derive  the  properties  of  the 
steady  state  and  the  equation  of  evolution  of  the  populations  of  the 
vibrational  levels,  where  we  look  at  the  effect  of  anharmonicity  by 
restricting  our  equations  to  a  two-level  system.  In  Section  VI  we  present 
some  numerical  results,  and  we  end  with  the  conclusions  in  Secion  VII. 

II.  Model 

The  internal  vibrational  mode  of  the  adsorbed  diatomic  molecule  is 
assumed  to  be  harmonic,  with  the  molecule  adsorbed  perpendicular  to  the 
surface.  The  infrared  laser  beam  is  collimated  and  polarized  in  such  a  way 
that  the  electric  field  is  parallel  to  the  axis  of  symmetry  of  the  molecule. 
This  simplification  is  valid  because  an  infrared  beam  can  only  create  an 
electric  field  perpendicular  to  the  surface  in  its  vicinity  [20].  The  laser 
is  tuned  in  resonance  with  the  internal  molecular  vibration.  The  adsorbate- 
substrate  electronic  system  is  taken  to  be  in  the  ground  state,  so  that  the 
vibration  of  the  adsorbed  molecule  creates  electron-hole  pairs.  This  causes 
a  strong  damping  of  the  vibrational  energy  into  the  metal.  This 
approximation  is  valid  if  ftu  »  kT,  where  T  is  the  temperature  of  the  metal, 
k  is  the  Boltzmann  constant,  and  u>  is  the  frequency  of  vibration  of  the 
molecule. 

The  Hamiltonian  can  be  written  in  the  second-quantization  form  as 

H  -  es  E*Eg  +  £[VaS  Eg  a  +  h.c.]  +  a+a  +  Vg(Oa+a  +  pbb+b 
S  S  b 

+  £  Vb(Ob+b  +  Hrv  +  Hr  +  KE  .  (1) 

b 
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Here  Eg  is  the  creation  operator  of  one  electron  in  the  state  |S>  of  the 

unperturbed  metal;  a^  and  are  creation  operators  of  electrons  in  the 

orbitals  |a>  and  |b>,  respectively,  of  the  free  diatomic  molecule;  £  is  the 

coordinate  associated  with  the  stretching  of  the  molecule;  and  KE  is  the 

kinetic  energy  associated  with  £.  |a>  is  assumed  to  be  the  only  molecular 

orbital  interacting  with  the  electronic  orbitals  |S>  of  the  metal,  and  this 

interaction  gives  the  factor  V  c.  e  +  V  (£),  e.  +  V.  (£)  and  e„  are  the 

ab  a  a  d  D  s 

Hartree-Fock  energies  of  the  electrons  in  the  orbitals  |a>,  |b>  and  |S>, 
respectively.  H  is  the  Hamiltonian  of  the  radiation,  given  by 

»r  ■  l  K  +  5>  ■  <2> 


Here  k  is  the  norm  of  the  wave  vector,  £,  of  the  photons,  where  k  is 
parallel  to  the  direction  of  propagation  of  the  beam.  It  »  kt^,  and  p£  is  the 
creation  operator  of  one  photon  of  wave  vector  k£g.  The  light  is  contained 

ITV 

within  a  cubic  quantization  volume  U.  H  is  the  interaction  Hamiltonian 
between  the  molecule  and  the  infrared  beam.  Since  we  are  only  concerned 
with  the  coupling  between  the  molecular  stretching  and  the  electronic 
system,  we  have  neglected  the  dependence  of  this  Hamiltonian  on  the 
positions  of  the  metal  nuclei  and  molecule  with  respect  to  the  surface.  We 
have  also  assumed  a  weak  dependence  of  V^g  on  versus  £.  This  is  a 
generalization  of  the  Hamiltonian  written  by  Persson  [10]. 

We  can  decompose  H  as 


H  -  H  +  H 


+  Hr  + 


HeV  + 


,rv 


(3) 


e 

where  H  is  the  electronic  Hamiltonian, 


H 


e 


l  eS  ESES  +  K s4‘  +  +  +  [  S  b+b  - 

S  S  b 


(4) 


and  can  be  formally  diagonalized  to  give 


H 


e 


l  «,  &  *  l  % ^ 

P  P 


(5) 


where 


a  “ 


£  <a|p> 
P 


(6) 


The  vibrational  Hamiltonian  is  given  by 

HV  =  <g|Va(C)a+a  +  £  Vb(C)b+b  |g>  +  KE  ,  (7) 

b 

where  |g>  is  the  ground  state  of  the  electronic  system.  Hv  is  assumed  to  be 
the  Hamiltonian  of  an  harmonic  oscillator. 


HV  «  -  #u>[V+V  +  |]  ,  (8) 

where  y  is  the  reduced  mass  of  the  molecule  and  V*  is  the  creation  operator 

rv 

for  a  vibrational  quantum.  H  is  then  given  by 

HrV  -  6CE  -  «((2^)1/2[V+  +  V])(-i(^)1/2  -  PkD  ,  (9) 

where  6  is  the  effective  charge  carried  by  the  stretching  coordinate  £, 

E  is  the  electric  field,  and  is  the  dielectric  constant  of  the  vacuum. 
Hev  is  given  by 

H®V  -  Va(U  a+a  +  £  Vb(Ob+b  -  <g|  y&U)  a+a  +  jvb(Ob+b|g>  .  (10) 

b  b 

We  approximate  by 

V  U)  -  TT  5  -  TT  <i>l/2<v+  +  v)  • 


(11) 


Ill .  Kinetic  Equation 

To  compute  the  evolution  of  the  vibrational  state  of  the  molecule,  we 
shall  use  the  Zwanzig  projector  technique  [19].  In  this  approach,  Zwanzig 
derived  an  exact  kinetic  equation  for  the  density  operator  of  one  subsystem 
interacting  with  another  one,  which  will  be  named  bath  subsystem.  In  our 
situation  the  first  subsystem  will  be  the  internal  vibrational  mode  of  the 
adsorbed  diatomic  molecule.  It  interacts  with  the  bath  subsystem,  which  is 
constituted  by  two  independent  subsystems:  the  laser  beam  and  the 
electronic  system.  The  interaction  Hamiltonian  between  the  vibration  and 
the  bath  is  given  by 

H1  -  C^)172  (V+  +  V)  B  »  (12> 

where 

B«Be®I  +Br$>I  .  (13) 

r  e 

Here  Ir  and  Ig  are,  respectively,  the  identity  operator  of  the  radiative  and 
electronic  subsystem,  §  symbolizes  the  outer  product,  and 

9e 

Be  -  •—  (aa  +  <g|aTa|g>)  (14) 


(15) 


GV 

We  have  removed  the  dependence  of  H  on  the  electronic  degrees  of  freedom 


corresponding  to  the  occupation  of  the  molecular  orbitals  |b>.  This 
simplification  comes  from  the  fact  that  these  orbitals  are  not  coupled  with 
the  electronic  orbitals  |S>  of  the  metal,  which  can  be  verified  by  computing 
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explicitly  the  memory  kernels  in  Eq.  (24)  or  (28)  below,  where  no 
specification  on  the  form  of  H*  is  needed. 

The  starting  operator  of  the  complete  system  [p(t  *  0)]  is  taken  to  be 
the  outer  product  of  the  vibrational  [pV(t  =  0)]  and  the  bath  (pb)  density 
operator. 


p(0)  =  pV(0)fcpb. 


(16) 


We  assume  the  vibrational  subsystem  to  be  in  its  ground  state  at  time  t  *  0, 


pV(0)  =  |0><0|  . 


(17) 


The  density  operator  of  the  bath  subsystem  is  itself  the  outer  product  of 
the  radiative  (pr)  and  electronic  (pe)  subsystem: 

Pb  »  P®  ®  Pr  .  (18) 

The  electronic  subsystem  is  taken  to  remain  in  its  ground  state, 

P6  =*  |g><g|  •  (19) 

The  density  operator  of  the  radiative  subsystem  corresponds  to  a  statistical 
superposition  of  coherent  states  with  different  phase  [21).  pr  is  then 
diagonal  in  the  radiative  Hamiltonian  eigenstate  representation, 


pr  *  \  Pkn  jknXnkj  . 
kn 


(20) 


We  note  that  p  (0),  pb,  HV  and  Hb  (=  He  Hr)  commute  with  each  other. 
Furthermore, 
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where  Tr  {  }  is  the  trace  over  all  the  electronic  and  radiative  states, 
er 

To  obtain  a  tractable  kinetic  equation,  we  shall  invoke  the  Born 
approximation  [12].  This  approximation  is  valid  if  the  interaction 
Hamiltonian  between  the  vibrational  and  bath  systems  is  small  in  comparison 
with  the  Hamiltonians  driving  their  evolutions  independently  (HV  »  H*  « 
Hb).  It  then  follows  from  that  the  kinetic  equation  can  be 
written  as  [12] 

i*pV(t)  =  [HV,  pV(t)]  -  |  ^  ds  (22) 

X  Trer^Hl*  exP("iHVs/tfMHI(-s),pV(t-s)&pb]exp(iHVs/rt)]}  » 

where  the  density  operator  of  the  vibrational  system  is  defined  by  Pv(t)  * 
Trer(p(t)},  and 

Hr(-s)  =  exp(-iHbs/yi)  H1  exp( iHbs/Hi) .  (23) 


Equation  (23)  can  be  written  in  terms  of  the  matrix  element  pV  (t) 

mn 


(*  <r'^v(t)|n>)  simply  by  multiplying  it  respectively  on  the  left  and  by  the 
right  with  <m|  and  |n>.  Thfc  derivation  in  the  case  where  m  *  n  is  given  in 
Ref.  12.  The  generalization  to  the  case  m  *  n  is  straightforward.  We 
observe  that  the  integrand  in  (22)  contains  a  double  commutator , 
corresponding  to  four  terms.  To  simplify  the  notation,  we  shall  track 
explicitly  just  one  of  these  terms,  keeping  in  mind  that  the  others  must  be 
treated  in  the  same  way.  We  then  arrive  at 


p!L(t)  *  -iu»  pv(t) 
nm  run  nm 


-  ~  f 


ds  (3  terms 


'  l  exp( 


iw  s) 
pn 


op 


where 


r-m  r,  rj  J 


r  ^  *•  r.  r  .  r. 


Hpm  *  <PlHl|n>  • 


<  K  >  =  Tr  {K-p  }, 
er 


a)  =  (EV  -  EV)/Ji  =  (n  -  m)u> 
run  n  m 


To  solve  these  coupled  equations,  we  transform  into  the  interaction 
picture  by  defining 


pY.(t)  =  P^  (t)  exp(-io)  t)  . 
run  run  run 


After  some  simple  algebra,  we  end  up  with 


p)L(t)  =  -  {3  terms  -  ]  exp[i(u  +  w  )t]  f  ds  exp(iu>  s) 

nm  wZ  £  po  nm  on 


x  <hL  (-s)>p'[  (t-s)}  .  (26) 

pm  no  op 

Without  any  interaction  (H^  =  0>),  the  quantities  pQp(t  '  s)  are  constant. 
With  H*  *  0^,  they  are  expected  to  vary  only  on  a  time  scale  much  greater 
than  the  oscillation  period  of  the  vibration.  Thus  the  only  terms  in  the 
sum  "^Qp"  which  give  rise  to  a  variation  of  t )  during  a  period 
corresponding  to  several  oscillations  are  the  ones  verifying 


w  +u>  =0  or  p  =  m  -  n  +  o  .  (27) 

po  nm  r 

Neglecting  all  the  others,  we  can  write  the  following  kinetic  equation  where 

« ““  <Hij(‘s)  ■&>  ■  <hL  HJj(-s)>*-  a  “L  - 

pL.(t)  =  -  [  ds 

nm  w2  J. 


t  J0 
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X  y{exp(iw  s)  <H*  H* _(-s)>  p^(t-s) 

£  no  no  on  nm 

o 

+  Iexp(i»mos)  <H^  hJb(-s)>  5^<t-s>]*) 
-^exp(l»0„s)  <H(m.n)+0(ln  5o,(»-n)+o(t-s) 


o 

+  texp(ia)  s)  <H 
om 


I 

(n-m)+o,n 


Hm  n('s)> 
in,  o 


(28) 


The  above  equation  gives  us  the  following  result:  only  the  elements  of 

the  vibrational  density  matrix  which  belong  to  the  same  diagonal  (p„  with 

|j  -  i|  =  constant)  are  coupled  in  the  kinetic  equation.  More  generally, 

even  for  an  highly  anharmonic  vibration  it  is  easy  to  see  that  the  diagonal 

elements  of  pV(t)  evolve  independently  from  the  others  [put  m  =  n  in  Eqs. 

(26),  (27)  and  (28)].  The  importance  of  this  result  will  be  emphasized  in 

Section  VII.  In  our  situation,  further  simplifications  can  be  made.  Using 

(12)  we  find  that  <H^  N,  H*  (-s)>  is  different  from  zero  only  if  o  ■ 

im-n;+o,m  no 

u  1/2 

nil.  In  the  case  where  o  ■  n  +  1,  it  is  equal  to  [(m+l)(n+l)] 
x  D(-s),  where  D(-s)  is  defined  as 

D(-s)  5  <B  B(-s)>  .  (29) 


The  kinetic  equation  can  then  be  written  as 

pY,(t)  * - - —  f  ds  { ( (n  +  1)  exp(-iojs)  D(-s)  p^(t-s) 

""  2#^  Jo  11111 

+  n  exp(iws)  D(-s)  p^Ct-s)) 

+  (m+l)[exp(-iujs)  D(-s)  pv  (t-s)]* 

nm 

+  m[exp(icos)  D(-s)pV  (t-s)]*} 
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-  {([(n+lXm+Ol^^expduis)  D(-s))  +  c.c.}  p^+1  m+1(t~s) 

-  {(vlrai  exp(-itjs)  D(-s))  +  c.c.}pV  ,  ,(t-s)  •  (30) 

n-i,m-i 

The  kinetic  equation  for  the  diagonal  elements  is 

P  (t)  *  -  -  Re  ds  {(n+1)  exp(-ius)  +  n  exp(ius)}  D(-s)  P  (t-s) 

n  ftp*  Jo  n 

-  (n+1)  exp  (iws)  D(-s)  Pn+^(t-s) 

-  n  exp(-iws)  D(-s)  Pn_j(t-s)  ,  (31) 

where  P  (t)  *  pV  (t)  *  pV  (t).  A  last  simplification  comes  from  the  fact 
n  nn  nn 

that 

<Be>  *  <Br>  =  5>  ,  (32) 

e  r 

where  <K>  =  Tr  {K*pe},  and  Tr  {  }  is  the  trace  over  all  the  electronic 

e  e  e 

states  (with  a  similar  definition  for  <  >r). 

The  memory  kernel  then  splits  into  two  parts,  one  related  to  the 
interaction  between  the  molecule  and  the  laser,  the  other  related  to  the 
interaction  between  the  molecule  and  the  substrate  [22], 

D(-s)  -  D*(-s)  +  Dr(-s)  ,  (33) 

where  De(-s)  »  <Be  Be(-s)>  ,  with  a  similar  definition  for  Dr(-s). 

e 

Because  at  time  t  ■  0  the  nondiagonal  elements  of  the  vibrational  density 
matrix  are  equal  to  zero,  they  will  remain  zero  at  all  times.  We  are  not 
interested  in  their  evolution  anyway,  because  the  mean  vibrational  energy 
does  not  depend  on  them. 


IV.  Computation  of  the  Memory  Kernel 


Let  us  first  compute  the  part  of  the  memory  kernel  due  to  HeV.  From 
(13),  (23)  and  (33),  we  obtain 

3e  -  .  . 

D  (-s)  =  (g^J  <g| (a  a  -  <g|aa|g>)exp(-iHes/lO 

x  (a+a  -  <gJa+ajg>)exp(iHes/J()|g>  .  (34) 


Using  (6)  and  its  Hermitian  conjugate,  we  can  transform  this  into 


De(-s)  =  (^)2  ^  ^  j<a|a>|2  |<a|p>|2exp[-(eQ  -  ep)s/K]  . 

ea  <  ef  eg  >  ef  (35) 

This  formula  has  been  derived  previously  by  Persson  [10].  The  new  feature 

which  we  shall  introduce  is  a  generalization  to  all  vibrational  transitions, 

the  justification  of  the  Markov  approximation,  and  the  formula  (39)  which  is 

not  clearly  established  in  [10]. 

Replacing  the  sum  by  an  integral. 


£  •*  J  de  p(c)  ,  p(c) j<aja>|2  -►  Pg(e)  , 


we  obtain 


De(-s)  «  (-gf)2  f  f  dea  J*  deg  pa(ea)  Pa(ep)  exp[-i(EfJ  -  co)s/*]  , 

0  Ef  (36) 

where  Pfl(e)  is  the  density  of  states  of  the  orbital  |a>.  If  we  assume  that 
the  width  of  pa(e)  is  on  the  order  of  leV  [23],  then  De(-s)  vanishes  for  s 
»  11/ leV  —  10  ^  s.  During  such  a  short  time,  the  system  cannot  evolve,  and 
D  (-s)  may  be  integrated  over  s  assuming  Pn(t-s)  z  Pn(t)  in  (31),  which 
is  the  Markov  approximation.  From  [10],  we  have 

Re  j  ds  exp(iws)  De(-s)  -  irrt2u>  • 
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Re  ds  exp(-iws)  D®(-s)  *  0  (37) 

J0 

Taking  into  account  the  fact  that  the  orbital  |a>  is  fourfold  degenerate 
[10],  we  end  up  with  a  rate  of  transition  from  the  vibrational  state  n  to  m 
equal  to 

R®  «7(nl  )o2  ,  (38) 

nm  4p  n,m+l  ' 

3c 

where  o  -  jf  pa(ef)  . 

As  shown  in  [10],  o  may  be  determined  experimentally.  It  is  the 

derivative  of  the  number  of  electrons  in  the  orbital  |a>  with  respect  to  the 

stretching  amplitude.  Persson  assumed  that  the  difference  between  the 

dipole  moment  of  the  adsorbed  (D  )  and  free  (D,)  molecule  is  due  to  this 

a  t 

charge  transfer.  Then,  the  effective  charge  (6)  associated  with  the 
stretching  of  the  adsorbed  molecule  is  equal  to  6  *  (6+  +  ido),  where  6+  is 
the  effective  charge  associated  with  the  stretching  of  the  free  molecule,  e 
is  the  charge  of  one  electron,  and  d  is  the  mean  distance  separating  the 
electron  in  the  orbital  |a>  from  the  metal  surface.  Using  the  fact  that  the 
transition  dipole  moment  of  the  free  and  adsorbed  molecule  are, 
respectively, 

D  .  s+(JL)l/2  „„  D  .8(JL)l/2  , 

f  2pui  a  2pu) 

we  can  estimate  o  to  be 


VDf  y2UU)\  1/2 
id  '  *  ; 


We  now  compute  the  part  of  the  memory  kernel  due  to  Hr  .  Some  simple 
algebra  leads  to 


Dr(-s)  "  2e^U  /  wk*<PkPk>eXp^‘i(1,kS^  +  ^<pk  pk>  +  1)exP(llJJks)}  •  (40) 

0  f* 


By  taking  U  large  enough,  we  may  replace  the  sum  by  an  integral 


<pk  pk>  n*“L*  * 


where  c  is  the  speed  of  light. 


We  also  neglect  the  term  "(+l)exp(-iaij[s),  which  gives  the  damping  of 
the  vibrational  level  due  to  the  vacuum  fluctuations  of  the  electric  field. 
Finally,  we  assume  the  spectral  intensity  of  the  laser  to  be  a  Lorentzian 
centered  around  u>. 


n(oL  )  *  N  - = - t—  ,  (41) 

ir[I”L  +  (u  -  ol^)  ] 

where  is  the  spectral  linewidth  of  the  laser  and  N  is  the  density  of 
photons.  Because  n(cj^)  vanishes  quickly  when  is  different  from  id,  we  may 
replace  in  (40)  by  id.  The  intensity  of  the  laser  is  given  by: 


0  nr,.x  BojN 
“  D  (0)  "  2/3 


Finally,  we  obtain  from  Eqs.  (40) -(42) 

r,[exp(iiD,  s)  +  exp(-ia»Ls)) 


$ 


[exp(ius)  +  exp(-iiDs)) 


<  'ft.t't.t j 


The  quantity  Dr(-s)  does  not  vanish  before  s  is  larger  than  1 . 
is  the  temporal  coherence  of  the  electric  field  and  may  be  larger  than  the 


characteristic  time  of  the  evolution  of  the  probabilities  Pn(t),  in  which 


case  the  Markov  approximation  is  not  applicable.  Furthermore,  we  may  not 
then  neglect  the  temporal  coherence  of  the  stretching  position  £(t).  To 


-1 


take  this  into  account,  we  replace  t  u  in  (31)  by  ±u»  +  if.  T  is  the 


temporal  coherence  of  £(t)  and  is  given  by 


10 


lift? 

4y 


(44) 


V.  Equation  of  Evolution  and  Property  of  the  Steady  State 
By  neglecting  the  integrals 


[  ds  exp(±i2u>s)  exp(-Ts)  P  (t- 
Jn  n 


s)  , 


we  can  express  the  kinetic  equation  (31)  in  the  simple  form 


*(t)  *  5>  ?(t)  +  2^ce0  J0  ds  exP('rs)  *(t's)  *  (45) 

where  r  *  rL  +  Tv,  ?(t)  is  the  vector  (Pg(t),  P^(t)  ...),  D*  is  the  matrix 


0  1 


0  -1 


0  0 


0  0 


-n  (n+1) 


and  E  is  given  by 


*  •  •  1  «  V*  »  "  »  *  -  *  •  ‘  «  *»"  *'*  * 


In  the  steady  state,  P(t)  is  zero  and  ?(t-s)  is  equal  to  ?(t),  such  that 
(45)  becomes 


?(t)  *  a[D>  +  pE*]  £(t)  *  (J 


(46) 


-ith  •  -  ^  and  »  -  ■ 

4|i  ,.2-2 

irn  uro  cGq 

In  fact,  by  eliminating  "■  3"  on  the  far  right-hand  side  of  Eq.  (46),  we 
obtain  the  kinetic  equation  within  the  Harkov  approximation.  In  the  steady 
state,  the  probabilities  P®  are  given  by 


Cl  ■  +  vl  Pn  -  K-l 


(47) 


This  kind  of  equation  admits  only  one  normalized  solution  [24],  given  by 


PS  -  (  -  P  -)n  — i- 
n  V1  +  p'  1  + 


(48) 


The  mean  number  of  vibrational  quanta  is  given  by 


<n>  «  jl  nP*  «  p  '  (49) 

n=0 


The  power  dissipated  into  the  metal  is  given  by  the  scalar  product  between 
aB>?s  and  (0,](oj, 21fao, . . .  ) ,  which  is  equal  to 


62I 

2pce0r 


Power  ■  |-aphw|  = 


molecule 


(50) 
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To  know  how  many  times  it  takes  to  the  system  to  reach  the  steady 
state,  we  must  integrate  Eq.  (45)  numerically.  This  leads  to  the 
iterative  equation 


P  (t+At)  *=  P(t)  +  At-a{[-nP  (t)  +  (n+1)  P  ..(t)] 
n  n  n  n+1 


+  pr 


[  ds  exp(-Ts)[nP 

J0  n 


^ ( t ”S )  +  (n+l)Pn+1(t-s) 


-  (2n+l)Pn(t-s)]}  . 


(51) 


To  iterate  this  equation  to  arbitrary  t,  we  need  not  know  all  the  history, 
Pn(t-s),  of  the  system  since, 

,t+AT 

j  ds  exp(-Ts)P  (t+At-s) 

J0  n 


-TAt 


-  e  (At'P  (t+At)  +  ds  exp(-Ts)P  (t-s)] 

n  n  1  n 


i: 


(52) 


When  the  anharmonicity  of  the  vibration  is  high,  the  laser  can  only 
couple  the  two  lowest .vibrational  levels.  The  interaction  between  the 
molecule  and  the  metal  entails  a  vibrational  transition  to  just  the  lower 
level.  Then  the  populations  of  the  levels  containing  more  than  one  quantum 
will  decay  rapidly  to  zero  (a  *  -  lifetime).  Furthermore,  if  at  t  *  0  these 
levels  are  not  occupied,  they  will  remain  empty  all  the  time.  In  that 
situation,  our  problem  reduces  to  a  two-level  system,  and  the  kinetic 
equation  reduces  to 

•  • 

?a(t)  ■  <x[(q  _j)?a(t)  +  pr(  _j)  j  ds  exp(-Ts)  ?a(t-s)J  .  (53) 

2 

Here  ctpr  «  0  /2,  where  n  is  the  Rabi  frequency  defined  by 

«2  “  fJ-|<l|M|0>|2/K2  . 

CC0 


where  M  =  6£  is  the  dipole  moment  operator,  and  ?a(t)  =  (Pa(t),Pa(t)) .  At 
the  steady  state,  the  population  of  the  upper  level  and  the  power  dissipated 
into  the  metal  are  given  by 


b  _ a 

2»  +  1  '  202  ♦ 


+  2or 


5ower  *  molecule  * 

Zp  T  1 


These  results  are  formally  equivalent  to  those  obtained  by  van  Smaalen, 
Amoldus  and  George  (see  Section  VI  of  Reference  16),  who  transformed  to 
dressed  states  to  compute  the  laser  heating  of  a  transparent  crystal  via 
adsorbed  atoms.  They  accounted  for  the  detuning.  A,  between  the  frequency 
of  the  laser  and  the  vibration,  and  they  neglected  the  lifetime  a  *  and  the 
laser  width  In  any  event,  by  taking  the  high  laser  intensity  limit,  we 

end  up  with  the  same  saturation  power 

PowerSat  =  molecule  *  (56) 

Comparing  (54)  and  (55)  to  (49)  and  (50),  we  observe  that  the  only  new 
feature  introduced  by  the  anharmonicity  is  the  saturation  factor  (2p  +  1) 
so  that  the  effect  of  the  anharmonicity  is  negligible  if  P  «  1. 

The  kinetic  equation  (53)  can  be  rewritten  in  terms  of  the  population 
inversion  as 


[F?(t)  -  Pa(t)]  =  -2aP 


a(t)  -  tJ2  J  ds  exp(-Ts)  [Pa(t-s)  -  Pa(t-s)] 


On  a  time  scale  shorter  than  T  ,  we  can  describe  analytically  the  behavior 

£ 

of  the  population  inversion.  Neglecting  the  term  -2oP^  and  estimating 
exp(-Ts)  as  unity,  we  can  write  the  solution  of  (57)  as 
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(P®(t)  -  pj(t)]  =  -  cos(nt)  .  (58) 

Here  we  have  taken  the  initial  conditions  to  be  Pq(O)  *  1  and  P®(0)  *  0. 

This  is  the  Rabi  oscillation,  and  such  behavior  cannot  be  revealed  using  the 
Markov  approximation.  This  oscillation  will  damp  on  a  time  scale  of  T  * 
to  reach  the  steady  state  defined  in  (54)  and  (55). 

VI.  Numerical  Results 

In  this  section  we  adopt  the  numerical  data  listed  in  Table  I,  which 

correspond  to  carbon  monoxide  molecule  adsorbed  on  a  copper  surface.  The 

laser  characteristics  correspond  to  an  infrared  beam  emanating  from  an 

optical  parametric  oscillator  pumped  by  a  YAG  laser.  This  beam  is  focused 

2 

to  obtain  a  section  of  1  mm  and  strikes  the  surface  with  a  small  grazing 
angle  of  incidence. 

From  (39)  and  (46),  we  obtain  (10] 

a  -  rv  «  6.4  x  10U  s'1  .  (59) 

-1  -12 

Since  the  temporal  coherence  of  the  position  £(t)  (Tv  *  1.5  *  10  s)  is 
much  smaller  than  the  temporal  coherence  of  the  laser  electric  field  (T^  *  > 
10  **  s),  we  may  make  the  estimation 

r  -  r  +  r.  -  r  -a  .  (60) 

v  L  v 

With  this  condition,  the  mean  number  of  vibrational  quanta  is  inversely 
proportional  to  the  square  of  a.  This  reflects  the  fact  that  this  number  is 
limited  by  the  rate  of  damping  of  the  vibrational  energy  into  the  metal,  as 
well  as  by  the  broadening  in  energy  of  the  vibrational  levels.  The 
broadening  decreases  the  absorption  efficiency  of  the  quasi  monochromatic 
laser  beam.  This  does  not  occur  if  the  linewidth  of  the  laser 
is  greater  than  the  width  of  the  vibrational  levels,  in  which  case 


Adopting  the  numerical  data  of  Table  I,  and  using  Eq.  (49),  we  obtain 
the  mean  number  of  vibrational  quanta  at  the  steady  state  to  be 

<n>  =  0  =  0.03  .  (62) 

The  power  transmitted  from  the  laser  beam  to  the  surface  is  calculated  from 
Eq.  (50)  to  be 

Power  »  7  *  10  ^  J  s  ^  molecule  ^  (63) 

From  (50)  and  (60),  we  observe  that  the  power  is  inversely  proportional  to 
a,  due  to  the  fact  that  Tt  is  much  smaller  than  T  .  If  the  interaction 
between  the  molecule  and  the  substrate  is  less  important,  we  have  a 
situation  where  TL  is  much  larger  than  in  which  case  the  power  is 
independent  of  a,  where  a  measures  the  strength  of  the  interaction  between 
the  molecule  and  the  substrate.  This  holds  if  the  hypothesis  of  harmonicity 
remains  valid,  which  means  that  0,  the  mean  number  of  vibrational  quanta  in 
the  steady  state,  is  not  too  high,  say  less  than  0.1.  The  region  where  our 
model  fails  depends  on  the  laser  characteristics  and  is  given  by 

log(f-)  <  log( - — )  ,  (64) 

L  *wrLMCC0 


where  we  have  assumed  a  <<  .  This  is  illustrated  in  Fig.  1. 

To  estimate  the  time  t  necessary  to  reach  this  steady  state,  we  compute 

the  evolution  of  the  populations  of  the  vibrational  levels  using  (51)  and 

(52).  We  adopt  the  initial  conditions  P  (0)  “6  The  evolution  of  the 

n  n ,  u 

mean  number  of  vibrational  quanta  and  the  occupation  probabilities  of  the 


three  lowest  levels  are  shown  in  Figs.  2  and  3  for  two  different  values  of 
(J.  The  populations  of  the  vibrational  levels  in  the  steady  state  for  ($ 
equal  to  1  are  presented  in  Fig.  4,  here  we  take  the  steady  state  to  be 
achieved  when  the  mean  number  of  vibrational  quanta  reaches  90Z  of  its  exact 
value  in  the  steady  state,  t  does  not  depend  noticeably  on  3  and  is 
estimated  to  be 


-1  -  -12 
t  —  4a  —  5. 6  *  10  s 


This  small  value  indicates  that  the  mean  number  of  vibrational  quanta  will 

follow  the  instantaneous  intensity  of  the  laser  since  the  intensity  varies 
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on  a  time  scale  of  the  order  of  10  s.  In  Fig.  5,  we  show  the  evolution  of 
the  mean  number  of  vibrational  quanta  both  with  and  without  the  Markov 
approximation.  We  observe  that  this  approximation  overestimates  the  rate  of 
evolution  of  the  system  and  gives  t  equal  to  2.3  o  * . 

VII.  Conclusions 

We  have  derived  a  simple  kinetic  equation  for  the  populations  of  the 
vibrational  levels  of  a  diatomic  molecule  adsorbed  on  a  metal.  The  molecule 
interacts  with  the  electronic  degrees  of  freedom  of  the  metal  and  an 
infrared  laser  beam.  This  equation  can  readily  be  solved  and  is  valid  for 
all  laser  linewidths,  infrared  absorption  widths  of  the  adsorbed  molecule, 
and  time  scales.  It  has  been  shown  that,  to  compute  the  evolution  of  the 
populations,  the  Markov  approximation  is  not  valid  if  the  temporal  coherence 
of  the  laser  field  is  greater  than  the  time  scale  during  which  the 
populations  are  changing.  This  failure  was  observed  earlier  by  Beri  and 
George  [22]. 


Our  kinetic  equation  corresponds  exactly  to  the  one  obtained  using  the 
random-phase  approximation.  The  use  of  this  simplification  has  entailed  the 
neglect  of  the  nondiagonal  elements  of  the  density  matrix  to  compute  its 
evolution,  and  this  approximation  has  remained  questionable  in  the  case  of 
the  interaction  between  a  laser  and  a  molecule  [12].  It  has  been  used 
extensively,  either  implicitly  in  application  of  the  golden  rule  [14],  or 
explicitly  [13,22].  Our  theoretical  development  proves  that  this 
approximation  leads  to  the  correct  kinetic  equation  for  the  populations  of 
the  vibrational  levels.  This  conclusion  should  emphasize  the  validity  of 
some  theoretical  predictions  already  given  in  the  literature. 

Our  results  also  show  the  equivalence  of  the  steady-state  properties 
computed  using  the  golden  rule  [14],  the  random-phase  approximation  [12,22], 
the  Markov  approximation  [13],  the  transformation  to  dressed-states 
technique  [15-17],  or  the  Zwanzig  projector  technique.  In  the  specific 
situation  considered  in  the  previous  section,  the  probability  of  finding  an 
adsorbed  molecule  in  its  first-excited  vibrational  level  is  0.03,  which 
suggests  that  resonant  excitation  of  a  molecule  adsorbed  on  a  metal  surface 
could  provide  a  means  of  laser-enhanced  surface  reactions. 


Acknowledgments 

This  research  was  supported  by  the  Fonds  National  de  la  Recherche 
Scientifique  of  Belgium,  Le  Ministere  Beige  de  la  Politique  Scientifique 
(S.P.P.S.,  Bruxelles),  the  Office  of  Naval  Research,  the  Air  Force  Office  of 
Scientific  Research  (AFSC),  United  States  Air  Force,  under  Contract  F49620- 
86-C-0009,  and  the  National  Science  Foundation  under  Grant  CHE-8519053.  The 


United  States  Government  is  authorized  to  reproduce  and  distribute  reprints 


for  governmental  purposes  notwithstanding  any  copyright  notation  hereon.  We 
want  to  thank  Dr.  H.  F.  Arnoldus  and  Dr.  S.  van  Smaalen  for  helpful 
discussions  and  for  a  critical  reading  of  the  manuscript.  AP  thanks  the 
State  University  of  New  York  at  Buffalo  for  its  hospitality  during  his  visit 
in  February-March  1986,  when  much  of  this  research  was  carried  out. 


*5*i 


a 


i 

•%} 

'S. 


References 

[1]  T.  J.  Chuang,  Surf.  Sci.  Rep.  3,  1  (1983). 

[2]  M.  Mashni  and  P.  Hess,  Chem.  Phys.  Lett.  77,  541  (1981). 

[3]  C.  T.  Lin,  T.  D.  Z.  Atvars  and  F.  B.  T.  Pessine,  J.  Appl.  Phys.  48, 

1720  (1977);  C.  T.  Lin  and  T.  D.  Z.  Atvars,  J.  Chem.  Phys.  68,  4233 
(1978). 

[4]  K.  S.  Gochelashvili,  N.  V.  Karlov,  A.  I.  Orchenkov,  A.  N.  Orlov,  R.  P. 
Petrov,  Y.  N.  Petrov  and  A.  M.  Prokhorov,  Sov.  Phys.  JRTP  43,  274 
(1976). 

[5]  H.  E.  Bass  and  J.  R.  Fanchi,  J.  Chem.  Phys.  64,  4417  (1976);  A.  V. 
Khmelev,  V.  V.  Apollonov,  V.  D.  Borman,  B.  I.  Nikolaev,  A.  A.  Sazykin, 
V.  I.  Troyan,  K.  N.  Frisov  and  B.  A.  Frolov,  Sov.  J.  Quantum  Electron. 
7,  1302  (1977). 

[6]  M.  E.  Umstead  and  M.  C.  Lin,  J.  Phys.  Chem.  82,  2047  (1978). 

[7]  T.  J.  Chuang,  in  Vibrations  at  Surfaces,  ed.  by  R.  Caudano,  J.-M. 

Gilles  and  A.  A.  Lucas  (Plenum  Press,  New  York,  1982),  pp.  573-577. 

[8]  B.  E.  Hayden  and  A.  M.  Bradshaw,  Surf.  Sci.  125,  787  (1983);  R.  G. 
Greenler,  K.  D.  Burch,  K.  Kretzschmar,  R.  Klauser,  A.  M.  Bradshaw  and 
B.  E.  Hayden,  Surf.  Sci.  152/153.  338  (1985). 

19]  H.-J.  Krebs  and  H.  Liith,  in  Proc.  Int.  Conf.  Vibrations  in  Adsorbed 


Lavers  (Jiilich,  1978):  KFA  Reports  -  Jiilich,  Conf.  26,  135  (1978). 

[10]  B.  N.  J.  Persson  and  M.  Persson,  Solid  State  Commun.  36,  175  (1980). 

[11]  T.  Rantala  and  A.  RosAn,  Phys.  Rev.  B  34,  837  (1986). 

[12]  S.  Efrima,  C.  Jedrzejek,  K.  F.  Freed,  E.  Hood  and  H.  Metiu,  J.  Chem. 
Phys.  79,  2436  (1983). 

[13]  J.  T.  Lin  and  T.  F.  George,  Surf.  Sci.  115,  569  (1982);  J.  T.  Lin,  X. 
Y.  Huang  and  T.  F.  George,  J.  Vac.  Sci.  Technol.  B  3,  1525  (1985). 

[14]  Z.  W.  Gortel,  H.  J.  Kreuzer,  P.  Piercy  and  R.  Teshima,  Phys.  Rev.  B 
27,  5066  (1983);  Z.  W.  Gortel,  P.  Piercy,  R.  Teshima  and  H.  J. 
Kreuzer,  Surf.  Sci.  165,  L12  (1986);  C.  Jedrzejek,  K.  F.  Freed,  S. 
Efrima  and  H.  Metiu,  Surf.  Sci.  109,  191  (1981);  J.  T.  Lin  and  T.  F. 
George,  Phys.  Rev.  B  28,  70  (1982);  B.  Fain  and  S.  H.  Lin,  Physica 
138B,  63  (1986). 


Kj 

[15] 

H. 

F.  Arnoldus,  S. 

[g 

in 

press . 

K 

[16] 

S. 

van  Smaalen,  H. 

3 

submitted. 

117] 

H. 

F.  Arnoldus  and 

AV  V  W  V  V  V  V  •  •  •  v' « * .  >  •  • .>  A 


[18]  C.  W.  Ford  and  W.  H.  Weber;  Surf.  Scl.  129,  123  (1983);  B.  N.  J. 
Persson  and  M.  Persson;  Surf.  Sci.  97,  609  (1980). 

[19]  R.  Zwanzlg,  J.  Chem.  Phys.  33,  1338  (1960);  R.  Zwanzig,  in  Lectures  In 
Theoretical  Physics.  Vol.  Ill,  ed.  by  W.  E.  Brittin,  B.  W.  Downs  and 
J.  Downs  (Interscience,  New  York,  1961),  p.  106  ff. 

[20]  R.  G.  Greenler,  J.  Chem.  Phys,  44,  310  (1966);  50,  1963  (1969). 

[21]  M.  Sargent  III,  M.  0.  Scully  and  W.  E.  Lamb,  Laser  Physics,  (Addison- 
Wesley,  Reading,  massachusetts,  1974),  Chapt.  XVIII. 

[22]  A.  C.  Beri  and  T.  F.  George,  J.  Chem.  Phys.  83,  2482  (1985);  J.  Vac. 
Sci.  Technol.  B  3,  1529  (1985);  Z.  Phys.  B  60,  73  (1985). 

[23]  For  an  example  of  CO  on  Ni(lll)  see  B.  Gumhalter,  Phys.  Rev.  B.  33, 
5245  (1986),  and  H.  J.  Freund  and  R.  P.  Messmer,  W.  Spiess  and  G. 
Wedler,  Phys.  Rev.  B,  33,  5228  (1986).  For  CO  on  copper  cluster,  see 
Ph.  Avouris,  P.  S.  Bagus  and  A.  R.  Rossi,  J.  Vac.  Sci.  Technol  B  3, 
1484  (1985). 

[24]  H.  F.  Arnoldus,  Compt.  Phys.  Commun.  33,  347  (1984). 

[25]  B.  N.  J.  Persson  and  R.  Rydberg,  Phys.  Rev.  B.  24,  6954  (1981). 


m 

:*v 

AW 


•»V»  .v 


N.  ■  .  .  4 


A 


/■/v 

•/’V 


r. 


•  W  * 


To  account  for  the  screening  of  the  electric  field  by  the  electronic 


polarizability  of  the  adsorbed  molecule,  we  adopt  6  =  0.6  e  [25]. 

Characteristics  of  the  Laser  Beam 

-9 

pulse  duration:  15  *  10  s 
-2 

pulse  energy:  10  J 


beam  cross  section:  10  ®  m^ 

surface  enhancement  of  the  intensity  by  constructive 
interference  between  the  incident  and  reflected  beam:  -  4  [20] 

12  -2  -1 

mea.  'ntensity:  -  3  «  10  J  m  s 
laser  width:  Acu^  <  1  cm  *  —  3  «  10^^  s  * 


Figure  Captions 


1.  Dependence  of  the  power  dissipated  into  the  metal  on  the  strength  (a)  of 
the  interaction  between  the  admolecule  and  the  substrate,  where 

62I 

Power  ■  lim  Power(a)  *  —  • 

■“  a-0  2rLMC£0 

The  region  where  the  hypothesis  of  harmonicity  fails  is  darkened  with 
points. 

2.  Evolution  of  the  occupation  probabilities  of  the  three  lowest 
vibrational  levels,  Pg(t),  Pj(t)  and  P^(t),  and  of  the  mean  number  of 
vibrational  quanta,  <n>,  for  3  ■  1. 

3.  Same  as  Fig.  2,  with  3  *  0.01. 

4.  Occupation  probabilities  of  the  vibrational  levels  at  the  steady  state 
for  3  ■  1. 

5.  Evolution  of  the  mean  number  of  vibrational  quanta  under  the  Markov 
approximation  (<n>^)  and  without  this  approximation  ^<n>nonj()  for  0  * 
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